Astronomy & Astrophysics manuscript no. radovich 


February 2, 2008 


(DOI: will be inserted by hand later) 





The VIRMOS deep imaging survey: III. ESO/WFI deep U-band 
imaging of the 0226-04 deep field. * 

M. Radovich 1 , M. Arnaboldi 1 - 2 , V. Ripepi 1 , M. Massarotti 1 , H. J. McCracken 3 - 7 , Y. Mellier 4 ' 5 , E. Bertin 4 - 5 , 

G. Zamorani 7 ' 8 , 

C. Adami 6 , S. Bardelli 7 , O. Le Fevre 6 , S. Foucaud 6 , B. Garilli 9 , R. Scaramella 10 , G. Vettolani 7 , A. Zanichelli 8 , 

E. Zucca 7 

o 

o . 

CNJ , 1 I.N.A.F.,Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131 Napoli, Italy 

• 2 I.N.A.F.,Osservatorio Astronomico di Pino Torinese, Via Osservatorio 20, 10025 Torino, Italy 

' Universita degli Studi di Bologna, Dipartimento di Astronomia, Via Ranzani 1, 40127 Bologna, Italy 

Institut d'Astrophysique de Paris, 98bis Bd. Arago, 75014 Paris, France 
Observatoire de Paris, LERMA, 61 Avenue de l'Observatoire, 75014 Paris, France 
6 Laboratoire d'Astrophysique de Marseille, Traverse du Siphon, 13376 Marseille Cedex 12, France 
I.N.A.F., Osservatorio Astronomico di Bologna, Via Ranzani 1, 40127 Bologna, Italy 
Istituto di Radioastronomia del CNR, Via Gobetti 101, 40129 Bologna, Italy 
Istituto di Astrofisica Spaziale e Fisica Cosmica del CNR, Via Bassini 15, 20133 Milano, Italy 
10 I.N.A.F., Osservatorio Astronomico di Roma, Via Osservatorio 2, 00040 Monteporzio Catone (Roma), Italy 

■>sj- . Received XXX; accepted XXX 

(N . 

Abstract. In this paper we describe the (/-band imaging of the F02 deep field, one of the fields in the VIRMOS Deep 
Imaging Survey. The observations were done at the ESO/MPG 2.2m telescope at La Silla (Chile) using the 8k x 8k Wide- 
Field Imager (WFI). The field is centered at a(J2000)=02 A 26"W and <5(J2000)=-04°30'00", the total covered area is 0.9 
deg 2 and the limiting magnitude (50% completeness) is Uab ~ 25.4 mag. Reduction steps, including astrometry, photometry 
and catalogue extraction, are first discussed. The achieved astrometric accuracy (RMS) is ~ 0.2" with reference to the /-band 
catalog and ~ 0.07" internally (estimated from overlapping sources in different exposures). The photometric accuracy including 
uncertainties from photometric calibration, is < 0.1 mag. Various tests are then performed as a quality assessment of the data. 
They include: (i) the color distribution of stars and galaxies in the field, done together with the BVRI data available from the 
VIMOS survey; (ii) the comparison with previous published results of {/-band magnitude-number counts of galaxies. 
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1 . Introduction The survey aims at covering ~ 16 deg 2 in UBVRI with a 5 

cr limiting magnitude Iab ~ 24.5, and a smaller 1.2 deg 2 area 

The VIMOS imaging survey (Paper I) represents the prepara- reaching ^ „ ^ . e ^ VIM0S deep fieW Tfl ^ v _ baadt 

tory step of the deep redshift survey, which is now being car- expected 5 a limking magnitudes are Uab „ 24 mag for the 

ried out with the VIMOS spectrograph at the VLT UT3 at wMe survey> ^ „ 26 mag for ±e deep suryey In ^ ^ 

Paranal, Chile, by the VIMOS consortium. This preparatory suryey contains Qyer 1q6 galaxies fa five colors and represents 

multi-wavelength imaging survey has been done at the Canada- a major ady£mce Qyer Qther preyiouf , deep multicolor surveySj 

France Hawaii Telescope (CFHT) for the BVRI bands (Le probing structures on scales of „ 2Qh -i Mpc at z ~ 1 . 



Fevre et al. 2003, Paper I hereafter; McCracken et al. 2003, 
Paper II hereafter), and at the ESO NTT telescope for the K 
imaging (Iovino et al. 2003, in prep). The near ultraviolet part 
of this survey was carried out at the ESO MPI 2.2m telescope 
with the WFI 8kx8k camera, in the framework of an approved 
ESO Large program (164.0-0089) scheduled since Period 63 in 
1999, both in visitor and service mode. 



Surveys in the radio continuum (Bondi et al. 2003) and with 
XMM (Pierre et al. 2003) were carried out on the deep field 
of the VIMOS imaging survey. The catalogs from the U and 
BVRI deep field are also foreseen to be used for the optical 
identification of the radio continuum survey (Ciliegi et al. in 
prep). 

Until the availability of data from ground-based (e.g. the 
Send offprint requests to: radovich@na.astro.it CFHLS and VST surveys) or space (e.g. the GALEX project) 

* Based on observations carried out at the ESO MPI 2.2m telescope surveys, wide-field (/-band data with medium depth and large 
located at La Silla, Chile. sky coverage are not available yet. The (/-band is essential 
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Fig. 1. The VIRMOS F02 deep field. The images of the three point- 
ings (PI, P2 and P3) are here combined for display purposes only: 
they have been treated separately in data reduction and catalog extrac- 
tion. The full field is displayed in the upper plot. The outer solid-line 
rectangle shows the position of the ~ 1.2 deg 2 field taken in the BVRI 
bands. Overlaid are the images of the three C-band pointings; masked 
areas at the borders of the images were removed. Note the circular 
shape of the PI pointing, due to vignetting in the Loiano filter. The 
region in the inner box, where PI and P3 overlap, is the one zoomed 
in the bottom plot: the accuracy of the astrometry can be seen from the 
absence of double sources. 



to study the star formation properties of field galaxies in the 
nearby universe, and is also particularly suited to identify star- 
burst galaxies, AGN, and QSO at moderate redshift. These 
aspects, as well the determination of the f/-band luminosity 
function and the analysis of the morphological properties of 
galaxies compared to other bands, will be the subject of sepa- 
rate papers. 

Compared to other existing U BVRI surveys, the VIRMOS 
F02 deep field offers a good compromise between covered area 
and depth. The ESO imaging survey (EIS) in the Chandra Deep 
Field South (Arnouts et al. 2001) is comparable in depth to 
our F02 deep field but covers a smaller area (~ 0.25 deg 2 ). 
The Canada-France Deep Field survey (McCracken et al. 2001) 
consists of four independent fields of ~ 0.25 deg 2 each, and 
the f/-band limiting magnitude is ~ 1 mag deeper than the 
VIRMOS deep field. The Combo-17 survey (Wolf et al. 2003) 



covers three fields with a total area ~ 0.78 deg 2 and is com- 
parable in depth to the F02 deep field in the t/-band. Deeper 
surveys (e.g. the William Herschel Deep Field, Metcalfe et al. 
2001 or more recently the FORS Deep Field, Heidt et al. 2003) 
were carried out on smaller areas (~ 7 x 7arcmin 2 ) only. 

In this paper we will describe the observations, photomet- 
ric calibration, catalogue extraction and validation for the deep 
field. The calibration steps are strongly related to those fol- 
lowed in the case of the overlapping deep BVRI filters, which 
are described in Paper II. We refer to this paper for a detailed 
discussion of calibration issues when they are the same and 
emphasize those aspects that are peculiar to the {/-band data. 

The paper is organized as follows. In Section 2 we describe 
the observations for the deep pointings, the observing strategy 
and data reduction issues (pre-reduction, astrometric and pho- 
tometric calibration). Catalog extraction is discussed in Section 
3. Validation of the catalogs through a number of tests (com- 
parison of stellar and galaxy colors with simulations, number 
counts) is done in Section 4. Conclusions are drawn in Section 
5. 

2. Observations and data reduction 

2.1. Observations 

The f/-band deep pointing was observed during several 
observing runs since 1999 with the aim of covering ~ 
1 deg 2 of the multiwavelength deep field planned within the 
VIMOS deep redshift survey. This deep field is centred at 
a(J2000)=02 /! 26 m 00 5 and (5(J2000)=-04°30'00"(see Fig. 1) 
and observations were carried out with the Wide-Field Imaging 
(WFI) mosaic camera mounted on the ESO MPI 2.2m tele- 
scope at La Silla, Chile. The camera was mounted at the 
Cassegrain focus of the telescope, giving a field of view of 
34x33 arcmin. It consists of a mosaic of 8 CCD detectors with 
narrow inter-chip gaps, yielding a filling factor of 95.9 % and 
a pixel size of 0.24". The WFI CCDs have a read-out noise of 
4.5 e~ pix~' and a gain of 2.2 e~ ADU 1 . 

The area covered for the U deep field consists of three 
pointings partially overlapping (see Fig.l). When this survey 
started in 1999, a standard {/-band Johnson filter was not avail- 
able at the ESO MPI 2.2m telescope, so our team borrowed 
the one available at the Loiano observatory, and used this till 
a suitable {/-band filter was acquired by ESO. One pointing 
was therefore imaged using the Loiano U filter, which is a cir- 
cular filter partially vignetting the field. Pointings 2 and 3 were 
imaged with the ESO U/360 filter. 

The filter transmission curves are shown in Fig. 2; central 
wavelengths and FWHMs are given in Table 1, together with 
covered area and exposure times for each pointing. Note that, 
compared to the ESO filter, the Loiano filter is shifted to the 
red and also has a red leak at A ~ 7000A. 

The efficiency during the observing campaign was ham- 
pered by bad weather conditions and a problem at the pri- 
mary mirror support system which caused strong astigmatism 
on some images. 

The exposure time for each individual frame was usually 
2000s, and 1000s when astigmatism of the ESO MPI 2.2m 
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Fig. 2. Transmission curves of the Loiano and ESO U/360 filters; 
curves before (dashed) and after (solid) the convolution with the CCD, 
telescope and atmosphere transmission are displayed. Note the red 
leak for the Loiano filter. 



telescope was severe. The total exposure time was 13.9 hrs for 
pointing PI, 20.3 hrs for pointing P2, and 17.8 hrs for pointing 
P3; see Table 1 for additional information on the observing log. 

For each pointing, the sets of exposures were acquired in 
a dithered elongated-rhombi pattern. This sequence of dithered 
exposures ensures the removal of CCD gaps in the final coad- 
ded image, better flat fielding correction, efficient removal of 
bad pixel and columns and a more accurate astrometric solution 
in the final coadded image. The total area covered by the three 
pointings is ~ 0.93 deg 2 . The effective area is actually smaller, 
~ 0.71 deg 2 , since: (a) a fraction of PI is overlapping with P2 
and P3; (b) the border regions of the images were masked when 
the noise was significantly higher due to either vignetting (PI) 
or the dithering steps; (c) a small fraction of the area in P2 and 
P3 is not covered in BVRI and therefore was not taken into 
account when catalogs were built. 

The atmospheric turbulence produced an average seeing of 
approximately FWHM ~ 1.3". 



2.2. Data reduction 

Pre-reductions were carried out using the MSCRED package in 
IRAF. The CCD mosaic frames were bias and dark corrected, 
and flatfielded. Flatfield images were constructed combining 
a whole series of twilight sky images, taken during each ob- 
serving night. Once the images were flatfielded, we noticed the 
presence of residual structures in the background sky, which 
should be treated in order to get a flat sky background in the fi- 
nal coadded images. A "super-flatfield" image was constructed 
from all the dataset for each pointing, using a 3-cr rejection al- 
gorithm, for the removal of sources in the field. The procedures 
for astrometry, photometry and coaddition were approximately 
the same as those adopted in the analysis of BVRI data. More 
details are given in Paper II; a short description follows where 
details peculiar to t/-band data are emphasized. 




2.42 2.44 
a 



Fig. 3. Astrometric residuals for unsaturated, point-like sources over 
the full field (PI, P2 and P3). Up: radial residuals; the inner and outer 
circles enclose 68% (cr = 0.21") and 90% (cr = 0.48") of the sources 
respectively; the cross shows the centroid, (<r ff , crj)=(-0.02", 0.03"). 
Bottom: residuals are displayed versus right ascension and declina- 
tion; the lines display the residuals per coordinate (68%: <r„ = 0.14", 
cr d = 0.16"; 90%: cr„ = cr 6 = 0.28"). 



2.2.1. Astrometry 

Astrometry was performed using the astrometrix tool, which 
is part of the WIFIX 1 package developed for the reduction of 
wide-field images. Astrometrix allows to compute an astro- 
metric solution using both an external astrometric catalog and 
the constraint that the position of overlapping sources in differ- 
ent CCDs must be the same (global astrometry). 

One of the main goals of this survey was to provide 
f/-band fluxes of the sources detected in the BVRI images. 
As described in more detail in Paper II, source detection for the 
BVRI images was done as follows. A single^ 2 image was first 
built from the BVRI images, source detection was then done 
using SExtractor (Bertin & Arnouts 1996) in dual mode. This 



1 The WIFIX package was developed by M.R. in coop- 
eration with the TERAPIX team; it is freely available at 
http://www.na.astro.it/~radovich. 
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Pointing 


Date 


Filter 


A c 


FWHM 


Area 


Ditherings 


Tot. exp 


Seeing 








(A) 


(A) 


(deg 2 ) 




(far) 


(arcsec) 


1 


Nov 99 


Loiano 


3620 3750 


527 460 


0.27 


25x2000s 


13.9 


1.4 


2 


Oct 00, Nov 00 


ESO U/360 


3404 3540 


732 536 


0.30 


22xl000s+ 25x2000s 


20.3 


1.3 


3 


Aug 01, Oct 01, Apr 02 


ESO U/360 


3404 3540 


732 536 


0.29 


32x2000s 


17.8 


1.2 



Table 1. Details concerning the pointings for the U deep field. Filter central wavelengths and FWHM are given before and after (in italics) the 
convolution with the telescope, CCD and atmosphere transmissions. The area includes only the region where sources have been included in the 
catalog: it does not include border regions with vignetting or a low signal to noise ratio. 




Fig. 4. Weight map produced after coaddition of the PI pointing. 

requires that source positions in all the images match at a sub- 
pixel level. Such level of accuracy was achieved as outlined in 
Paper II. The astrometric solution was first computed for the 
/-band image taking the USNO A-2 (Monet et al. 1998) as the 
astrometric reference catalog. A catalog of sources was then 
extracted from the resampled /-band image and used as ref- 
erence catalog for the other bands (UBVR). During the astro- 
metric procedure, the offset of {/-band detected sources with 
respect to those matched in the /-band catalog was first com- 
puted for each CCD separately. This allowed to correct the dis- 
placement introduced by atmospheric refraction from the /- 
to the U- band. Finally, in the global astrometry step the as- 
trometric solution was constrained for each CCD by both the 
positions from the /-band catalog and those from overlapping 
sources in all the other CCDs. Even if the three pointings were 
then resampled and coadded separately, in the global astrome- 
try we used positions from all of them to increase the astromet- 
ric internal accuracy by using the position of the same source in 
the overlapping region of two different pointings (see Fig. 1). 

To compute an accurate astrometry for i/-band data is 
somewhat more difficult than for BVRI data because of the 
smaller number of bright stars in the t/-band. This implies that 
the astrometric solution may heavily rely on extended sources, 
for which the emission may peak at different positions in differ- 
ent bands. The achieved RMS astrometric accuracy, measured 
using /-band selected point-like sources only, is < 0.2" (see 
Fig. 3). The internal RMS, computed from overlapping sources 
in different exposures, is much smaller, ~ 0.07". 



2.2.2. Photometry 

Photometric calibration was done taking several Landolt fields, 
so that we had at least one Landolt star in each CCD. Landolt 
fields were imaged also in B,V to derive the filter color term, 
for each f/-band filter used. 

The calibration equation becomes: 

U-u = a(U-B)+B-eX, (1) 

Landolt fields were taken at different airmasses (X) dur- 
ing the night for all pointings with the exception of P2; (U-B) 
colors of the Landolt stars used in the calibration are in the 
range < (U-B) < 1.5. This allowed us to compute the 
extinction coefficient e: a value close to that expected for La 
Silla was obtained. For P2 this was not possible, so we took 
the average La Silla extinction coefficient for that period. Zero 
points, color terms and extinction coefficients computed for 
each pointing are given in Table 2. Recent photometric investi- 
gations (Manfroid, Selman & Jones 2001) showed that the zero 
point changes across the WFI mosaic with a maximum differ- 
ence of ~ 0.08 mag as a consequence of non-uniform illumi- 
nation and scattered light. As the Landolt fields are not wide 
enough to cover the whole WFI area and hence to compute one 
zero point for each CCD independently, the zero points we ob- 
tained represent an average over the whole mosaic from the 
different observed Landolt fields. The photometric accuracy in 
our data is therefore limited to ~ 0.05 mag, corresponding to 
about half the maximum difference in zero points over the en- 
tire WFI mosaic. 

From the knowledge of the system transmission curves we 
then computed the AB corrections, that are displayed in the last 
column of Table 2. 

To account for changes in airmass and atmospheric trans- 
parency, we used the Photometrix package in WIFIX. This 
package first applies the astrometric solution found for each 
CCD and then looks for overlapping sources in different CCDs. 
It then computes/br each mosaic (defined as the set of 8 CCDs 
taken for each exposure) an additional term to the zero point 
(z r ,i) such that the average differences in fluxes of overlapping 
sources are minimized. An exposure taken in photometric con- 
ditions was chosen as reference, so that for it z r j=0: this im- 
age was taken in the same night as the Landolt fields used to 
make the photometry calibration. Since the P2 and P3 point- 
ings are not connected, and PI was taken with a different filter, 
we needed to run this step for each pointing separately. 
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Pointing 


Date 


Filter 


a 


P 


£ 


RMS 


X 


AB corr 


1 


Nov 1999 


Loiano 


0.21±0.02 


22.15 ± 0.06 


0.43±0.02 


0.06 


1.125 


0.73 


2 


Nov 2000 


ESO U/360 


0.08±0.02 


22.01 ± 0.02 


0.49 


0.06 


1.105 


0.94 


3 


Oct 2001 


ESO U/360 


0.08±0.02 


22.10 ±0.02 


0.44±0.03 


0.04 


1.139 


0.94 



Table 2. Zero points, color terms and airmasses for the different pointings and filters. Date and airmass (X) refer to the image taken as reference 
when the relative photometry was computed. The last column gives the values to be added to the Vega-system magnitudes to obtain AB 
magnitudes. 



2.2.3. Coaddition 

After the astrometric solutions and flux scaling factors were 
computed, image coaddition was done using the SWarp tool 
developed by E. Bertin 1 . SWarp allows to subtract the back- 
ground, resample the images according to the astrometric solu- 
tion, apply the flux scaling factor, and finally combine them. As 
in the case of the BVRI images, resampling was done using a 
"Lanczos-3" interpolation kernel, which corresponds to a sine 
function multiplied by a windowing function. The coaddition 
was done computing the median of the images to optimize the 
rejection of spurious sources as cosmic rays or satellite tracks. 
As a consequence of the use of two different filters for the PI 
and P2, P3 pointings, it was not possible to produce a single 
coadded image for the whole field. We therefore produced one 
image for each pointing; the size, pixel scale (0.205" pix~') and 
orientation of each image was the same as that of BVRI data so 
that catalogs can be extracted with the same x 2 technique (see 
Paper II). This allows to extract catalogs using SExtractor in 
dual mode, with the x 2 image as reference. In SWarp, each 
image may be associated to a weight-map to properly weight 
pixels during the coaddition step. This was particularly useful 
in the case of the Loiano filter, due to vignetting in the outer 
regions. Weight maps were first created using normalized flat 
fields; pixels flagged in the Bad Pixel Maps provided for WFI 
by ESO 3 were then set to in the weight map. For the point- 
ing PI, we also set to those pixels with a value < 0.6 in the 
normalized flat field, thus allowing us to remove in the coaddi- 
tion the regions affected by vignetting. Fig. 4 shows the coadded 
weight map for PI. 

3. Preparation of the catalogs 

3.1. Photometric properties 

A first estimate of the detection limiting magnitude for point- 
like sources was done on the background RMS map: this map 
was extracted for each pointing using SExtractor and the me- 
dian background RMS (<x) was then computed. The limiting 
magnitude at 3<x and 5<x levels is mag(ncr) = B-2.5 \og(ncr y/A) 
where n = 3,5, A is the area of an aperture whose radius is 
the average FWHM of point-like sources, B is the zero point 
(Eq. 1). We obtained U A b ~ 25.8 at 5cr and Uab ~ 26.4 at 3cr. 

Because of their definition, these values for the limiting 
magnitudes refer to the unrealistic case of objects with a flat ra- 

2 SWarp is part of the TERAPIX software suite, available at 
http://terapix.iap . f r/sof t 

3 http://www.eso.org/science/eis/eis_soft/softJndex.html 




22 24 26 28 



U AB 

Fig. 5. Ratio of detected to input sources as a function of input magni- 
tude: the input catalog consists of simulated point-like sources with 
a flat magnitude distribution. The dotted lines show the adopted 
50% and 90% completeness magnitudes for each pointing. The his- 
tograms in the insets show the magnitude distribution of the spurious 
sources detected in the images with background only and no simulated 
sources. 

dial surface brightness profile. A more accurate determination 
of the photometric limits of the survey was then done as fol- 
lows. Catalogs were first extracted using SExtractor, where 
source detection was optimized using simulations of point-like 
sources. A background image was computed as outlined in 
Arnaboldi et al. (2002); a population of 2000 point-like ob- 
jects with the same PSF of bright isolated stars in the field was 
added on this image, following a flat magnitude distribution 
(20.5 < mag < 27) and a random spatial distribution: the tasks 
in the ARTDATA package in IRAF were used. SExtractor was 
then run on that image and it was measured how many objects 
were detected and the faintest magnitude reached for a given 
detection threshold. The number of connected pixels required 
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Pointing 


^comp 


"Jlim 


/*comp 






(90%) 


(50%) 


(90%) 


(50%) 


1 


25.0 


25.5 


25.7 


26.2 


2 


24.9 


25.3 


25.8 


26.3 


3 


24.9 


25.3 


25.7 


26.2 



Table 3. Completeness and limiting magnitudes and surface bright- 
nesses (mag arcsecr 2 ) derived from simulated point-like and extended 
sources respectively. All values are in the AB system. 



19 30 31 33 
U360 



19 30 31 
U360 



for source detection was set to 9 pixels: this gives a minimum 
signal to noise ratio of ~ 3<x for a detection threshold ~ 1 . The 
same process was repeated 20 times, giving a total of 40 000 
sources and around 1500 sources per 0.25 magnitude bin. At 
the same time, we monitored the number of spurious sources, 
by matching the output catalogue from SExtractor with the 
input catalogue. The optimal threshold was the one that mini- 
mized the number of spurious detections with Uab < 26, with- 
out loosing input sources. Following these tests, we set the de- 
tection threshold to ~ 0.7cr for the PI pointing and ~ 0.9cr for 
P2 and P3, over at least nine connected pixels. 

Fig. 5 shows the results of the simulations (detected/input 
sources). In Table 3 we list the magnitude of completeness (> 
90% of the modeled objects in the input catalog are retrieved) 
and the limiting magnitude (50% of the modeled objects re- 
trieved). The depth of PI is comparable to that of P2 and P3 
even if the exposure time is lower, since the Loiano filter is 
redder than the ESO filter and is thus somewhat more efficient, 
after convolution with the atmosphere. In order to check how 
many spurious detections should be expected at different mag- 
nitudes, we then run SExtractor on the background image be- 
fore simulated sources were added. The result is displayed in 
the histograms in Fig. 5, which show that there are no spurious 
detections for magnitudes brighter than the 90% completeness 
magnitude; the peak occurs for magnitudes fainter than the de- 
tection limit (Uab ^ 26) 4 . Residual spurious detections for 
magnitudes close to the 50% limiting magnitude are later re- 
moved by the cross correlation with the catalog derived from 
the^ 2 image. 

Similar simulations were done using extended sources, 
to compute the completeness and limiting surface brightness. 
Galaxies with De Vaucouleurs and exponential surface bright- 
ness laws were used as input sources to simulate elliptical and 
spiral galaxies respectively; total magnitudes of these sources 
were in the range 20 < Uab < 26. The same procedure as above 
was then followed, but this time we measured the peak surface 
brightness (defined as the surface brightness of the brightest 
pixel) of the simulated and detected sources. The so obtained 
values are also given in Table. 3. 



4 Spurious detections with magnitudes in the range 28 < Uab < 30 
are mainly due to correlated noise. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

U360-B U360-B 



Fig. 6. Difference in aperture magnitudes (5"diameter) of point-like 
and extended bright (18 < U < 22) sources taken with both the ESO 
(U360) and Loiano (UL) filters. 6-band fluxes were taken from the 
VIRMOS catalog. The upper and middle panels show the differences 
in magnitudes before (left) and after (right) the color correction is ap- 
plied (squares: P2; dots: P3). In the bottom panel (U - B) colors from 
the Loiano and ESO filters are compared. 



3.2. Catalog extraction 

In order to easily combine the catalogs derived from the U and 
BVRI images, it was decided to use the BVRI x 2 image cata- 
logue as reference for the extraction of sources in the {/-band. 

The BVRI x 2 image catalogue is significantly deeper than 
the catalogue which can be derived using only the f/-band im- 
age; the relative depths of the two catalogues can be judged by 
comparing the surface densities of objects, ~ 10 5 per square de- 
gree in the^ 2 image catalogue (Paper II), and ~ 10 4 per square 
degree in the U band catalogue. On the other hand, we would 
measure positive (spurious) flux in the background-subtracted 
U frame at the position of sources detected in the x 2 image 
even if no significant U source is detected (by definition a pos- 
itive flux would be measured in about 50% random positions). 
To remove these spurious detections we adopted the following 
strategy. For each pointing (PI, P2 and P3), we extracted two 
catalogs: 
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- One catalog was extracted using SExtractor in the "dual 
image mode", taking as input catalogue the one obtained 
from the^ 2 image of the BVRI dataset (see Paper II). 

- One catalog was extracted in single-image mode, using a 
filtering Gaussian kernel and the detection threshold set as 
discussed above. 

For both catalogs, Kron (magauto) magnitudes and aperture 
magnitudes at diameters of 15 pixels (~ 3") and 25 pixels 
(~ 5") were computed. As it concerns the required number of 
contiguous pixels, we used the same value that was adopted for 
the BVRI images (see Paper II). 

A cross-correlation was then done between coordinates of 
sources in the two catalogs; a maximum distance of 2" was 
adopted for matching. In addition, aperture magnitudes (5" di- 
ameter) were also used as a further constraint in the cross- 
correlation: matching sources were rejected when the differ- 
ence of //-band magnitudes obtained from the x 2 an d single 
image catalogues was > 0.5 mag. 

After the cross correlation we computed the average differ- 
ence between the magnitudes measured in the single and x 2 
images: for all the three pointings we obtain (AU) < 0.02 
mag with RMS values which change from < 0.05 mag for 
bright sources (U AB < 22) to ~ 0.1 mag for fainter magni- 
tudes. Considering (1) the different depths of the {/-band and 
X 1 images and (2) the fact that for extended sources the peak of 
the emission may be different in the i/-band compared to the 
BVRI bands, we decided to keep in the final catalog the magni- 
tudes measured on the single image rather than those computed 
using the parameters from the^- 2 image. 

As a consequence of the dithering strategy, regions close to 
the border of the image are covered by a small number of expo- 
sures: spurious sources are therefore produced either by a bad 
rejection of cosmic rays or by the high noise. We therefore re- 
moved from the catalogs all the sources located within a given 
distance from the border. //-band magnitudes were assigned 
to sources in the x 1 catalog only when they were detected in 
the single-image catalog with a magnitude brighter than the 
3cr limiting magnitude, U A b = 26.4. For magnitudes brighter 
than the completeness limit, unmatched sources are mainly 
false detections (e.g. halos around bright stars) or sources not 
found in the ^catalog (e.g. because of cosmetic defects). For 
fainter magnitudes, spurious detections due to correlated noise 
are more important in the single-image catalog (see Fig. 5) but 
are removed by the cross correlation with the x 2 catalog. The 
requirement on the magnitude helps to reject these spurious de- 
tections, removing < 6% of sources for U < 25 mag, ~ 10% 
for fainter magnitudes. 

The separation of extended vs. point-like sources was done 
on the /-band catalog, and is described in Paper II: such clas- 
sification was done in the range 18 < Iab < 21, that is between 
the /-band saturation limit and the magnitude beyond which 
the separation between resolved and unresolved sources is less 
reliable. We verified on the half-light radius versus //-band 
magnitude plots that sources classified as point-like fall on the 
expected locus. The same plots were also used to estimate the 
average seeing in each pointing that is given in Table 1 . In addi- 




-0.5 0.5 1 1.5 



Fig. 7. Comparison of colors for bright ( Uab < 23.5) point-like 
sources (dots) with: the library from Pickles (1998) (stars); the 
Kurucz model atmospheres (Kurucz 1979) with logg = 4.0, [M/H]=0 
(lower curve) and [M/H]=-5.0 (upper curve). Magnitudes in PI were 
trasformed to the ESO filter photometric system. 

tion, bright nearby galaxies which are saturated in the /-band 
but not in the //-band were flagged as extended. 

Finally, we computed for each filter the correction to bring 
the magnitudes to the AB system (see Table 2) and the Galactic 
extinction correction. The F02 field is characterized by a low 
interstellar extinction (see Paper I); according to the Schlegel 
et al. (1998) maps 5 0.023 < E(B - V) < 0.040, with an average 
value E(B - V) ~ 0.027 ± 0.002 mag. This translates to an 
average correction for extinction A v ~ 0.14 mag. 

3.3. Color correction of the Loiano filter to the ESO 
filter 

After catalogs were extracted, for each pointing we assessed 
the problem of bringing the //-band photometry to a single 
photometric system. We decided to convert the magnitudes 
measured in PI with the Loiano filter to the ESO photomet- 
ric system, since the latter is closer to the Johnson system. The 
advantage of doing so, rather than using Eq. 1 to transform both 
magnitudes to the Johnson system and then computing the in- 
ternal offset, is given by the fact that we can later use the ESO 
system transmission curve e.g. for comparison with model pre- 
dictions. 

As shown in Fig. 1, the P2 and P3 fields taken with the ESO 
filter partially overlap with PI. For the overlapping sources, 
we therefore have photometry for both the //-band filters; in 
addition, B-band photometry is available from the CFHT data. 
We first extracted from the catalog bright (20 < U < 22) point- 
like sources; 35 and 46 sources were found for P2 vs. PI and P3 
vs. PI respectively. From these we solved Eq.l with no airmass 
(e = 0), taking the ESO filter as reference (magnitudes not 
including the AB correction were used). Similar results were 
obtained for both P2 and P3: hence we merged the two sub- 
catalogs and finally obtained a = 0.178, B = -0.08 in Eq.l. 

5 The dust map was downloaded from: 
http://astron.berkeley.edu/davis/dust/data/data.html. 
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( B - V )ab 

Fig. 8. Comparison of simulated and observed (B— V) colors vs. / mag- 
nitudes. Top: color-magnitude diagram. The upper panel shows the 
synthetic colors expected in the F02 area. The lower panel shows the 
observed colors for {/-band selected (Uab < 23.5) point-like sources 
from the three pointings; the cuts at the faint and bright end are due to 
the fact that the selection of point-like vs. extended sources was done 
in the range 18 < Iab < 21 only. The solid line indicates the separation 
between halo and disk/thick disk stars. Bottom: The same as above, 
but plotted as a histogram. 



The color term, as expected, is intermediate between the values 
found for the two filters from the calibration with standard stars 
(see Table 2). In Fig. 6 we show the differences in magnitudes of 
both point-like and extended sources in the overlapping areas 
of the two filters before and after color correction. The average 
offset in magnitude after the correction is -0.01 mag (PI vs. P2) 
and 0.04 mag (PI vs. P3), with an RMS -0.1. The same plot 
also provides an indirect check that photometry in P2 and P3 is 
consistent, as each of them is consistent with the PI photometry 
in the overlapping areas. 

The consistency with magnitudes obtained in the Johnson 
system by applying the color terms listed in Table 2 was 
checked as follows. Both the ESO and Loiano magnitudes were 
first transformed to the Johnson system. The same procedure 
as above was then followed to compute the offset between the 
Loiano and ESO photometry: we found a = as expected and 
B = -0.1 mag. Loiano magnitudes were then transformed to 
the ESO system as discussed above, and then to the Johnson 




0.5 1 1.5 0.5 1 1.5 



(U- B )ab ( U - B )ab 

Fig. 9. (U — B) AB vs. (V — 1)ab colors for extended sources with 
Uab < 25. Colors are displayed as a grey-scale map where the level 
is proportional to the number of sources in a given color bin. Model 
tracks for early to late type galaxies are also displayed (starting from 
right to the left); the positions of lower and upper redshift are indi- 
cated by the filled and open circles respectively. Magnitudes in PI 
were transformed to the ESO filter photometric system. The four pan- 
els show sources in different bins of I AB , as follows, (a) I AB < 18, 
model tracks: 0.0 < z < 0.2. (b) 18 < Iab < 20, model tracks: 
0.0 < z < 0.4. (c) 20 < I AB < 22, model tracks: 0.0 < z < 0.8. 
(d) 22 < Iab < 24, model tracks: 0.0 < z < 1.5. The value in the 
top left corner is the fraction of galaxies detected in the (/-band with 
respect to those detected in the /-band, in each bin. 

system. The agreement is within ~ 0.08 mag in the range 
-2<(U-B)<2. 

4. Data quality assessment 

In order to check the data quality of the [/-band photometry, 
we performed a series of tests where we (i) compared colors of 
point-like and extended sources to those expected for stars and 
field galaxies; (ii) checked the number counts of galaxies and 
compared them with existing data in literature. In this anal- 
ysis, magnitudes in the PI pointing were first transformed to 
the ESO filter photometric system as described in Sec. 3.3; for 
those sources which were observed with both the Loiano and 
the ESO filter, the latter magnitude data were taken. In the final 
catalog, magnitudes are corrected for Galactic extinction. 

4.1. Stellar colors 

We first compared the colors of bright (Uab < 23.5) point- 
like sources with those obtained from the convolution of stellar 
templates (Pickles 1998) with the system transmission curves 
(filter + telescope + CCD + atmosphere). No correction for 
Galactic extinction is applied in this case. The result is dis- 
played in Fig. 7: it can be seen that the agreement is good 
only for (B — V) > 0.7. This is expected since for bluer colors 
the low-metallicity stellar population from the halo dominates 
over the disk population with solar metallicity. We therefore 
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computed the colors produced using the Kurucz model atmo- 
spheres (Kurucz 1979) with logg = 4.0 and [M/H]=0.0, -5.0 
to reproduce solar and sub-solar metallicities (see also Lenz et 
al. 1998). As displayed again in Fig. 7, the apparent excess in 
(U - B) is in agreement with the colors produced by the low- 
metallicity model. The absence of observed stars at redder col- 
ors ((U — B)> 2.5) is due to the cut in the (/-band magnitudes, 
as it is explained later in more detail. Most of the sources out- 
side the stellar locus with an ultraviolet excess (((/ - B) < 0.6) 
are likely to be quasars at z < 2.2. 

In order to verify whether the color distribution we obtained 
is consistent with that expected in our Galaxy, we show in Fig. 8 
(top) the synthetic color-magnitude (CM) diagram ((B - V) vs. 
/) for a region of sky centered on PI and a size of 0.9 deg 2 
(Degl'Innocenti & Cignoni, Private Communication). The code 
and recipes described in Castellani et al. (2002) were used to 
this aim: the code provides BVRI magnitudes for a Galactic 
stellar population model, including halo, disk and thick disk 
components. 

We selected from the simulations those sources with 18 < 
Iab < 21 as this is the range where the point-like classifica- 
tion was done in our data. We expect to find only halo stars 
(i.e. metal poor objects) around (B - V) ~ 0.5, a mix of halo 
and thick disk stars in the interval 0.5 < (B - V) < 1 .4, and 
mainly disk stars (metal rich objects) around (B - V) ~ 1.4. 
Fig.7 shows that (U - B) > 2 when (B - V) > 1: therefore, 
when we compare simulated and observed color distributions 
we need to take into account the selection introduced by the 
cut in the {/-band magnitude, Uab < 23.5. As the models do 
not provide (/-band magnitudes, we proceed as follows. The 
limits (U - B) > 2 and U AB < 23.5 imply that B AB < 21.5. 
From (B - V) > 1 we obtain V AB < 20.5. Finally, Fig. 15 in 
Paper II shows that (V - I) > 1 when (B -V) > 1 : we therefore 
obtain Iab < 19.5. 

The same diagram is plotted in Fig. 8 for the point-like 
sources in the VIRMOS deep field which were detected in the 
i/-band. An exact match with the number of disk stars can not 
be expected as the selection due to the cut in the {/-band was 
taken into account in an approximate way. The position of the 
two clumps centered at (B - V) ~ 0.5 and (B - V) ~ 1.4 is 
clearly visible (Fig 8, bottom). This confirms that the (U - B) 
excess seen in Fig. 7 is due to the halo stars in the F02 field. 

4.2. Field galaxy colors 

(U - B)ab and (V - I)ab colors were compared with those ob- 
tained from model galactic spectra convolved with the system 
transmission curves. This was done using the make_catalog 
tool which is part of the HYPERZ package (Bolzonella et al. 
2000): we used the built-in set of templates for early to late- 
type galaxies computed from the GISSEL98 spectral evolution 
library of Bruzual & Chariot (1993). As in Paper II, we divided 
our catalog of (/-band detected extended sources in different 
bins of 7 A b magnitude. Late-type galaxies are expected to in- 
creasingly dominate in the color distribution as we go to fainter 
magnitudes: early-type galaxies with faint I magnitude are too 
faint to be included in the {/-band catalog. 
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Fig. 10. Galaxy magnitude-number counts from the VIRMOS 
(/-band deep field compared with existing literature data, (/-band 
magnitudes are in Vega system. The filled circles show the VIRMOS 
data (statistical error bars are displayed but are too small to be seen); 
magnitudes from PI were transformed to the ESO filter photomet- 
ric system. Literature data are taken from the N. Metcalfe's col- 
lection. The dotted line shows the best fit computed in the range 
18.5 < U < 22.5. The dot-dashed lines show the 90% and 50% com- 
pleteness limits. The histogram displays the number counts corrected 
for incompleteness (see text for details). 

(i) Iab < 18. The brightest galaxies are clearly dominated by a 
low redshift population, z < 0.2, of early to late-type galaxies 
(Fig. 9a). 

(ii) 18 < I A b < 20. In Paper II it is shown that (B - V) vs. 
(V - I) colors for these galaxies are in agreement with those 
expected for a population of galaxies with z = 0.0 - 0.5. The 
median redshift expected for this magnitude range is z ~ 0.2, 
according to spectroscopic surveys like the CFRS (Crampton 
et al. 1995). Consistent results are found looking at the (U - B) 
vs. (V - 1) colors (Fig. 9b), where early to late type galaxies are 
all seen to contribute to the observed colors. 

(iii) 20 < Iab < 22. The colors are well reproduced by a popu- 
lation dominated by late type galaxies with z ~ 0.0 - 0.8, with 
a peak at z ~ 0.35 (Fig. 9c). 

(iv) 22 < I A b < 24. The limiting magnitude Uab ~25 implies 
that in this range we mainly select galaxies with bluer colors 
((U - B)ab < 0.2): most of these are likely to be late type 
galaxies at z < 1, with some contribution from higher redshift 
galaxies (Fig. 9d). Note that in this magnitude range less than 
30% of the galaxies are detected in the f/-band. 

4.3. Number counts 

The comparison of number counts of galaxies with literature 
data provides a good check for both the efficiency of the star- 
galaxy separation (done in the /-band in our case) and the qual- 
ity of our photometry. For (/-band data, this is complicated by 
the few published number counts and by the heterogeneity of 
the data (e.g. different photometric systems, no correction for 
galactic reddening). A collection of (/-band number counts 
is provided in the Vega-system by N. Metcalfe's (http://star- 
www.dur.ac.uk/~nm/, see also Metcalfe et al. 2001); we se- 
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lected only those measurements from CCD observations. They 
are displayed in Fig. 10. The offset among the different mea- 
surements is due to both the different {/-band filters and the 
absence of reddening correction for some of the literature num- 
ber counts. For example, in the case of Metcalfe's WHDF data 
E(B - V)whdf ~ 0.02, so that after dereddening the WHDF 
number counts should be shifted by — 0.1 mag (Heidt et al. 
2003). 

The f/-band VIRMOS number counts were normalized 
to the area covered by the unmasked regions, taking into ac- 
count the overlapping regions (~ 0.7 deg 2 ). Number counts not 
corrected for incompleteness are displayed as filled circles in 
Fig. 10. Our data are in very good agreement with those of 
Arnouts et al. (2001), which were also taken with ESO/WFI. 
A least-squares fit in the range 18.5 < U < 22.5 gives a 
slope d(log N)/dm = 0.54±0.06. The histogram in Fig. 10 shows 
the number counts corrected for the ratios of input to detected 
sources displayed in Fig. 5 (Sec. 3.1): after this correction, 
number counts are in agreement with the extrapolation from 
brighter magnitudes up to U ~ 24.5 (U A b ~ 25.5). 

5. Summary 

f/-band data obtained in the framework of the VIRMOS 
preparatory imaging survey for the F02 deep field were pre- 
sented in this paper, as a complement to the BVRI data pre- 
sented in Paper II. Observations, data reduction and catalog ex- 
traction issues were first discussed. Various quality assessment 
tests were then performed. A good agreement is found between 
the observed stellar colors and those computed from the Kurucz 
model atmospheres; a component with sub-solar metallicity is 
required to fit the halo stellar population which dominates over 
the disk population for (U - B) < 1.5. This is confirmed by the 
comparison with the color distribution of stars in our Galaxy 
based on the Castellani et al. (2002) model. The colors of ex- 
tended sources were compared with those obtained from tem- 
plate spectra of galaxies. These tests also allowed to check the 
photometric consistency of f/-band and BVRI data. Number 
counts for extended sources were finally computed and com- 
pared with other f/-band data in literature. 
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